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Abstract: Using holographic QCD based on D4-branes and D8-anti-D8-branes, we 
have computed couphngs of glueballs to hght mesons. We describe glueball decay 
by explicitly calculating its decay widths and branching ratios. Interestingly, while 
glueballs remain less well understood both theoretically and experimentally, our re- 
sults are found to be consistent with the experimental data for the scalar glueball 
candidate /o(1500). More generally, holographic QCD predicts that decay of any 
glueball to At^q is surpressed, and that mixing of the lightest glueball with qq mesons 
is small. 
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1. Introduction 

Glueballs, as excitations of gauge-invariant composite operators in Yang-Mills theo- 
ries, remain illusive. Although the existence of the glueballs (of various types, such 
as scalar glueballs and tensor glueballs) is expected, their experimental identification 
in the hadron spectra remains difficult.* This difficulty is largely due to the inability 
to compute reliably couplings of glueballs to ordinary mesons in strongly coupled 
QCD. Lattice QCD predicts for the mass of the lightest scalar glueball to be around 
1600-1700 MeV but it doesn't yet provide information on the glueball cou- 

plings and decay products/widths, which are indispensable for their identification. 
The Large Hadron Collider (LHC) will likely yield a huge amount of hadronic data, 
which can lead to progress in revealing the mystery of glueballs. 

In this paper, we explicitly compute the couphngs between light glueballs and 
light qq mesons, by using holographic QCD. AdS/CFT (gauge/gravity) correspon- 
dence (duality) P, ^ is one of the most important developments in string theory, 
and, holographic QCD refers to the application of AdS/CFT to QCD studies. The 
basic claim of the AdS/CFT correspondence is that correlation functions of gauge- 
invariant composite operators in large Nc gauge theories at strong 't Hooft coupling 

*For details, we refer readers to the section of "Non-qq candidates" in Meson Particle Listings 
in and to Q for recent discussions. 
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correspond to classical gravitational computations in higher dimensional gravity the- 
ories in curved backgrounds. The correspondence has been applied to (i) computation 
of glueball spectrum in large Nc pure Yang-Mills theory and to (ii) qq meson spec- 
tra/dynamics in large Nc QCD, which we review briefly below. These efforts have 
been quite successful in reproducing lattice and experimental data of hadrons, even 
though the real QCD is recovered in the "CFT" side only when one incorporates 
various corrections in the large Nc and large 't Hooft coupling expansion. Here we 
combine these two efforts, (i) and (ii), in order to calculate couplings between the 
glueballs and the qq mesons, in the large Nc QCD. 

The key merit of using holographic QCD is the fact that one not only can calcu- 
late the hadron spectra, but also can compute explicitly their couplings. It provides 
a more powerful method for constraining these couplings than the chiral perturba- 
tion technique. In particular, since glueballs are expected to be heavier than 1 GeV, 
derivative expansion in chiral perturbation becomes unreliable. Furthermore, current 
lattice calculations are not well suited for computing dynamical quantities such as 
decays and couplings. Our paper represents a first principle calculation for glueball 
decays, though in the approximation where the holographic duality is valid. 

Let us briefly review here the holographic study of the two sectors (i) and (ii). 
Glueball studies began at the early stage of AdS/CFT correspondence, since they 
should exist in pure Yang-Mills theories whose supersymmetric version was the ba- 
sic building block of the correspondence. Witten was the first to suggest a reliable 
way in breaking the supersymmetries thus allowing one to treat the bosonic Yang- 
Mills theory 0. The gravity dual is the near horizon limit of a classical solution 
of 10 dimensional type IIA supergravity representing A^^^ D4-branes wrapping an 
5*^ with anti-periodic boundary condition for fermions. After various developments 
along this direction p, |T0|, |1T|, a complete spectrum of scalar/tensor glueballs in 
four-dimensional Yang-Mills theory was given in |T^ where they appeared as gravi- 
ton/dilaton/tensor fluctuations in the Witten's gravity background. The calculated 
glueball spectrum is consistent with the lattice computations, and in this paper, we 
compute glueball decays, based on this spectrum [jl2| . The lightest glueball is a scalar 
state with quantum numbers J^*" = 0+"*". Since this state should be the easiest one 
to identify in hadronic data, we concentrate here on this lightest scalar glueball for 
explicit computations. The gravity dual of this lightest scalar glueball corresponds 
to a specific combination of metric fluctuations. 

The quark sector, (ii), is obtained in AdS/CFT correspondence by the intro- 
duction of flavor D-branes []T3| intersecting with color Nc D4-branes. Various D- 
brane configurations (and also phenomenological holographic models) describing fla- 
vor/chiral physics have been proposed (see 0, |15|, |16|, [l^, |l^] for a partial list). 



Among these, we shall use the Sakai-Sugimoto model |Tj, |T8[, which has been quite 
successful in reproducing various facets of low energy QCD dynamics while maintain- 
ing its string-theoretical origin. The Sakai-Sugimoto model uses Nf D8-branes and 
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Nf anti-D8-branes as the flavor D-branes, and their intersection with color Nc D4- 
branes gives rise to string excitations corresponding to the quarks. Among various 
ways of introducing flavor D-branes, the Sakai-Sugimoto model beautifully realizes 
spontaneous chiral symmetry breaking and chiral dynamics in QCD.^ The resulting 
low energy theory for the quark sector is dual to the probe D8-brane worldvolume the- 
ory (higher dimensional Yang-Mills theory) in the Witten's supergravity background. 
In particular, the qq mesons are described as Kaluza-Klein (KK) decomposed mass- 
less fields on the probe D8-brane. 

We would hke to compute the couplings between the glueballs and the qq mesons 
in this setting. In the dual description through the AdS/CFT, they correspond to 
the supergravity fluctuations and the Yang-Mills fluctuations on the D8-branes, re- 
spectively. These two sectors are coupled in the combined system of supergravity 
plus D8-branes. We substitute the fluctuations (wave functions) of the supergravity 
fields (corresponding to the glueball) and the D8-brane massless flelds (mesons) into 
the D8-brane action and integrate over the extra dimensions, to obtain the desired 
couplings. Combining sectors (i) and (ii) is important not only due to its phenomeno- 
logical impact but also because this represents the first computation in holographic 
QCD of the couphngs between the supergravity fiuctuations and the fields on the 
probe D-branes. 

Once the couplings are obtained, we can compute the decay widths for various 
decay channels of a glueball, and study its possible mixings. Because the whole qq 
meson sector is combined into the D8-brane action, which is a higher-dimensional 
Yang-Mills lagrangian, several interesting mesonic features follow.-^ For example, at 
the leading order (in the expansion of the large 't Hooft coupling), glueball decay to 
Attq is prohibited. There is no direct Attq coupling to the glueball, and, furthermore, 
glueball - p meson coupling also does not allow the Atiq decay mode. As for the 
mixing, we can show that the lightest glueball has no mixing with qq mesons at 
the leading order of our expansion. These are our main predictions based on the 
holographic QCD. 

The organization of this paper is as follows. In section ||, we explicitly compute 
the glueball couplings in the holographic QCD, and obtain the interaction lagrangian. 
We study generic features of the glueball decay following from the holographic QCD. 
In section ^, we compute the decay widths based on these interactions. We list 
possible decay products, and obtain widths for various allowed decay modes. We 
next compare these with the experimental data. It has been argued that /o(1500) is 
the most plausible candidate for the lightest scalar glueball and we find that our 
results are consistent with the hadronic data for /o(1500). We reproduce the narrow 
width of the /o(1500), and also the decay products/branching ratio, qualitatively. 

'''Although the quarks of the model are massless (and so the pions are massless), it does not 
present a problem for our purpose. For recent discussion on obtaining massive pions, see p9| . 
^One is the reproduction of the vector-meson dominance as shown in uM. 
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In section we provide a summary, discussions, and a list of directions for future 
studies. 

2. Glueball Interaction 

Holographic QCD, in particular the Sakai-Sugimoto model, has provided a novel 
unified view of the mesons in a multi-flavored QCD. All the mesons appear just 
as KK decomposed massless fields living in higher dimensions. As a consequence, 
there exist many interesting relations among couplings between the mesons. The 
Skyrm term is one example. Here our concern is with the glueballs, which live 
in a different sector in the dual side, i.e. in the supergravity fluctuations, not on 
the flavor D-branes. But the holographic features found in the meson sector are 
inherited also for the glueball-meson couplings. This is because these couplings are 
also controlled by the flavor D-brane action, and shares the same flavor structure 
(commutator structure of the non-Abelian massless fields on the D8-branes). In this 
section, after reviewing the dual descriptions of the glueballs and qq mesons, we 
describe the generic features for glueball decays dictated by the holography. Finally, 
using the holographic QCD, we explicitly derive the couplings between the lightest 
scalar glueball and light qq mesons. 

2.1 Brief review of holographic QCD: glueballs and mesons 
2.1.1 Glueball sector 

Glueballs are gauge-invariant composite states in Yang-Mills theory, and their duals 
are fluctuations in near horizon geometry of black-brane solutions. There are nu- 
merous ways to break supersymmetries by deforming the AdS^ x solution with 
which the original AdS/CFT correspondence was derived. Among them, the Wit- 
ten's background |^ is "reliable" in the sense that it knows how the supersymmetries 
are broken in the Yang- Mills side: anti-periodic boundary condition for the fermions 
on the D4-brane worldvolume. 

Let us review briefly the description of the gravity dual for the lightest glueball 
in the four-dimensional QCD. It corresponds to supergravity fluctuations in the Wit- 
ten's classical background in 10 dimensions. The type HA supergravity solutions can 
be written conveniently in the notation of the 11 dimensional supergravity, and in 
that notation the Witten's solution is a doubly Wick-rotated AdS-j blackhole, 

ds" = ^ {f{r)dT^ + Vf.udx^dx") + ^{f{r))-'dr'' + h^dnj (2.1) 

where /(r) = 1 — R^/r^. L and R are the parameters of the solution, /i, i/ run 
from to 4, and the direction is the 11th dimension (the M-theory circle). The 
r direction is compactified to a circle, and its radius is fixed as L^/ (3i?) so that the 
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background is non-singular, and the manifold is smooth, around the "end" of the 
spacetime solution (at r = R). 

The S'^ part is not necessary in the following discussion, so we integrate that 
part (and also the flux) to obtain the M-theory supergravity action reduced to 7 
dimensions: 

S = ^j^V, j d'x^-deiGMN [R{Gmn) + ] . (2.2) 
Here V4 = 87r^/3 is the volume of a unit 5*^, and we have followed the notation of 



0.§ 

In ||12[, a complete bosonic spectrum was given, and the lightest state has quan- 
tum numbers J^^ = O"*"^ in terms of the ■ • ■ spacetime. The metric fluc- 
tuations for this lightest state in the action ( ^.2[ ) were explicitly obtained in [|T^^" 
as 

V - ^ {',.^,H^r) - (i + ^^,H(r)) (2.4) 

= M2L2(10r6-4i?6)2 ^^''^^^^^"^^' = -j2fH{r)G{x). 

Here (7(0;°, ■ ■ ■ ,a;^) is the glueball field in the real 1+3 dimensional spacetime, and 
M is the mass of the glueball.** 

The mass squared for this glueball state was found in [|T^] to be = 7.308i?^/ L"^, 
by solving the eigen-equation following from the equation of motion of the 7-dimensional 
AdS supergravity ( ^.2] ) given by the mass-shell condition of the glueball field (□ — 
M^)G = 0, 

For our later purpose, it is useful to change the coordinate to a dimensionless Z 
defined by 

r/R = K^/^, K = l + Z^. (2.6) 



5 Further integration of the r and the x'^ (M-theory circle) gives a 5 dimensional AdS gravity 
action, 



S = ^^fnlJ,L2u j drd^x r^V/V-detCMTV (^R{Gmn) + ^ ] (2.3) 



30 

Here M, N run through (0, 1, 2, 3, r), and we have used the M-theory ^ type IIA relations Ru — gsh 
and 2kI, = {2^fllgl 

^For the analogue state in three-dimensional QCD, see 

"in there is a typo in equations (38) and (40) (the sign of the functions c and h is opposite) 
**The excitation tower for these graviton-dilaton fluctuations is denoted as ^4 in [l2|. 
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Z = corresponds to the bottom of the background r = R, and the branch Z{> 0) 
is smoothly connected to the branch Z{< 0) [l^. In this Z coordinate, the eigen- 
equation becomes 

The appropriate boundary condition for solving this is 

= 0, H{Z = 0)^0, H{Z = oo) = 0. (2.8) 



dZ^ 



z=o 



Other fluctuations (such as the state corresponding to 2++ glueball) can be con- 
structed in the same manner |T^ . 

2.1.2 qq meson sector 

The dual of the quark sector is the probe flavor D-branes intersecting with the color 
D-branes. To make sure that we are not simply constructing a phenomenological 
model but "deriving" the hadronic interactions from the first principle, we need to 
follow so-called top-down approach from string theory. The Sakai-Sugimoto model 



1^ is the best known top-down construction for multi-flavor quarks. 

The qq mesons are described in the Sakai-Sugimoto model by the flavor D8-brane 
action in the Witten's background (written with type IIA string metric), 

Sm = -(27ra')'rD8Tr J d'x e-y^detg ^ f-g'^^FpgFRs + ^chcm- Simons- (2.9) 

Here gpq is the metric induced on the D8-brane worldvolume spanning the directions 
0, 1, 2, 3, r, S*^, and we have already expanded the Dirac-Born-Infeld (DBI)action to 
the second order in the Yang-Mills field strength.* The normalization of the genera- 
tors of the gauge group is chosen as Tr T^Tfe = 6ab- The parameters in the background 
metric ( p.l|) are related to the notation of the Sakai-Sugimoto model by L = 2_Rss5 
R = 2\/ RssUkk and ( p.6| ), where Rss denotes "i?" in the original papers of Sakai 
and Sugimoto [|1^, 0. The typical mass scale appearing in all the computations in 

O, [g is Mkk = mpl/KRss^'. 

Again, we integrate out the irrelevant 5*^ part, leading to 



(2.10) 



* There is a tadpole for closed string modes, but in this paper we neglect the back-reaction to 
the metric due to the presence of the D8-brane because it will be small in the large Nc limit. 
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where z = UkkZ and /i, z/ = 0, 1, 2, 3. (The Chern-Simons term in ( p.9|) will be irrel- 
evant to our discussion of the glueball decay; see section ||.) The KK decomposition 
along z (equivalently r, Z) in [jl^, |l8l is 



A, = Mz)^ixn, = M^)p>.i^n (2.11) 

where we suppress all the other higher components, since we are interested in the 
decay of the glueball to light mesons. The eigenfunction ipi should satisfy the eigen- 
equation which follows from the action ( p.lOj ), 

-K'^^dziKdz^i) = \iipi (2.12) 

where Ai = 0.669 for normalizable V'l; leading to the mass squared for the p meson, 
rn^p = XiMj^j^. The pion is massless, with its eigenfunction (po{z) ocl/K. 

The trace in (|2.1CI|) is for the matrix- valued mesons. For Nf flavors, the pions and 
the p mesons are Nf x Nf matrices. The overall trace part of the pion should yield 



a mass from the chiral anomaly (see [T^, for the description of the supergravity 



counterpart) but it is not included here. In this paper we take Nf = 2, hence the 
overall trace part of the pion is rj (or 7]') meson, while that of the p meson is u. 
Except when adjoint indices are explicitly written, we include r] [r]') and u in the 
matrix notation of the fields vr and p. 

2.2 Generic features of holographic glueball decay 

Our strategy to compute the interaction between the glueballs and the qq mesons is 
very simple. Since we know how all these hadrons are described in the dual side (as 
in (p^ ) and (|2.11|) ), we substitute them into the D8-brane action ( |2.9| ) and integrate 
it over the extra dimensions. In the original Sakai-Sugimoto model, the induced 
metric g in ( ^.9| ) was just the background metric, but now the glueball appears as a 
fluctuation in the induced metric and the dilaton in the D8-brane action. 

Since the appearance of the glueball doesn't break the non-Abelian structure 
of the D8-brane action, we can expect that some generic features of the glueball 
coupling may be read in the D8-brane action. As an obvious check, the glueballs 
should be flavor-blind; This can be seen as the fact that the supergravity flelds are 
gauge invariant with respect to the gauge transformation on the D8-brane. As a 
consequence, couplings of the glueballs to qq mesons are universal against flavors. 

Gauge invariance in higher dimensions also constrains the meson interactions. 
For example, as Sakai and Sugimoto have shown, the Skyrm term in the pion self- 
interactions is encoded in the structure of the higher dimensional Yang-Mills la- 
grangian. In our case, even though glueballs (=gravity and dilaton fluctuations) 
are now included, this flavor structure is almost unaltered. We note the following 
interesting features: 

(a) There are no glueball interactions involving more than two pions. 
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(b) For glueball coupling to p and pion, the p meson couples to the pion as if the 
pion were charged under the p meson gauge field. 

(c) Direct coupling of a glueball, G, with more than five mesons are suppressed by 
large 't Hooft coupling. 

(a) is easily seen by noting the fact that tt appears in but there are no {A^)"^ terms, 
with n > 2 in the D8-brane action ( |2.9| ). (b) follows from the fact that ~ [A^, Az\ 
can be decomposed as [p^, c?^7r]. These two features are precisely what were observed 
in the Sakai-Sugimoto model for the pure qq meson sector, and now inherited to 
the glueball couplings. Finally, (c) is due to the fact that the higher dimensional 
Yang- Mills action (|2.9|) does not have A^ term. For (a) and (c), DBI corrections give 
terms but they are suppressed by a' or equivalently the large 't Hooft couphng. 
To be more explicit, we can list the couplings which appear in ( p.9|) : 

GTr(7r2), GTr(7r, [vr, p]), GTr([7r, p]^), GTr(p^), GTr(p[p,p]), GTr([p,p]2). (2.13) 

Here we omit the derivatives and also possible indices. There are in fact no other 
couplings, and this is a generic result from the holographic QCD for interactions 
involving a single glueball. Even for multi-glueball vertices, this fiavor structure is 
maintained. 

A direct consequence of this fiavor structure is the fact that glueballs cannot 
decay to Atiq. The decay channel to Atiq appears in the F^^ term which is in the 
higher DBI corrections and thus suppressed by the large 't Hooft couphng. So, the 
holographic QCD predicts that, among the decay products of the glueballs, Atiq is 
suppressed. 

It is important to note that we work here in a "holographic gauge" (|2.11|) of 
the D8-brane action, in which interactions are seen in the simplest and the most 
transparent way. On the other hand, in the = gauge Jl^ the broken chiral 
symmetry is manifest because the pions fields appear in the action as [/ = exp(z7r//^). 
Thus this gauge is appropriate for comparison with the chiral perturbation theory, 
but the interactions are complicated. A different gauge choice in the D8-brane action 
leads to a different field definition of the four dimensional fields and, of course, this 
does not change the physics. Therefore, there should be some hidden structure in the 
QCD effective action, at least in the large N^. limit, because of the higher dimensional 
gauge symmetry. For example, the vector meson dominance is its consequence [0]. 
The suppression of the Attq in the glueball decay is also a consequence of this hidden 
structure which is not manifest in a gauge choice other than the "holographic gauge" 

For explicit computations, we concentrate on the couplings of the lightest scalar 
glueball, because of its phenomenological interest. But it is obvious that couphngs 
of other glueball excitations can be computed in the same manner. One of the 
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phenomenologically interesting excitations is the 2++ state. Fortunately, the super- 
gravity fluctuation for this 2++ state is simple (see for example whose notation 



we follow), in particular it consists of only the fluctuation of the metric components 
of /i, z/ = 0,1,2,3. So it in fact couples to the 4 dimensional part of the energy- 
momentum tensor of the D8-brane Yang-Mills action. Therefore the coupling should 
be of the form 



(2.14) 



at the quadratic order in vr and p^. Here F^i, is the field strength of the p meson at 
its linear order, F^;, = dfj,p^ - d^p^. 

Also the second lightest 0"''+ glueball is simple, since the supergravity fluctuation 
is involved with the dilaton that has a very simple coupling to the D8-brane Yang- 
Mills fields. We expect that this kind of simplicity in the glueball couplings may 
give some constraint on the decay products and decay widths, in particular the spin 
dependence of the decay product, and may serve as a smoking gun for identifying 
the higher glueball states in the meson spectroscopy.^ In the next subsection, we 
explicitly compute the interaction lagrangian of the lightest scalar glueball. 

2.3 Interaction of the lightest scalar glueball 

In this subsection we derive the interaction lagrangian of glueballs with light qq 
mesons (the pions and the p mesons). First, we need to fix the normalization of 
eigenfunctions in higher dimensions, in the dual side. Then, we substitute all the 
fluctuations into the D8-brane action and perform integration over extra dimensions 
to obtain the interaction lagrangian. The interaction lagrangian includes a possible 
mixing between the glueball states and the qq mesons; however, we show that there 
is no mixing for the lightest glueball. 

2.3.1 Normalization of the fluctuation fields 



Sakai and Sugimoto 17, 18 have determined the normalization of the eigenfunctions 



ipi{z) and (po{z) for the p meson and the pion respectively: 

\R%'ul/irM9:\2nar j dZ K-''\i^,f = 1, (2.15) 
!^ss'f^KK^D8^4^7;' (2W)2 j dZ K{<f)o)' = 1. (2.16) 



2 
3 
3 



With these, substituting ( |2.11| ) into the D8-brane action ( p.9| ) (the metric is fixed 



with its background value), we obtain canonically normalized kinetic terms for the 



''^In fact, we will find in the next subsection that there is no mixing between the lightest 0+"'" 
glueball (or the lightest 2++) with qq mesons, at the leading l/^/Wl order. 
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p meson and the pion, 

^0 = -Tr j d'x {^(9,vr)2 + \fI^ + ^AiMi^pjj . (2.17) 

The field strength F is that of the p meson, F^^ = d^pi, — d^p^. 

On the other hand, the normahzation of the fiuctuation eigenfunction H{Z) for 
the glueball has not been carried out in the past. The normahzation of H{Z) in 
( p.4|) should be fixed in such a way that substitution of the expressions ( p. 41) into the 
supergravity action (p.2|) results in, after the integration of the extra dimensions, a 
canonical kinetic term for the glueball field G{x), 

S = const. - j (^-{d^Gf + ^M^G^^ + 0{G^). (2.18) 

We numerically solve the differential equation ( |2.7| ) for H{Z), to get-^- 

j d'x^/-det Gmn (^RiG) + 

= -0.0574^ {H{Z = 0)f J d^xdrdx^ [d^.Gd^'G + M^G^] . (2.19) 

So, using the expressions for the M-theory gravity coupling /tn, and also the relations 
between the supergravity parameters and the QCD parameters [|17|, 

^^-T^' t'KK ^ H^^WJ. (2.20) 
we obtain the normalization 

{H{Z = 0))"^ = 0.00978gYMN^^^ Mkk. (2.21) 
2.3.2 Glueball interaction lagrangian 

Once the normalization of the eigenfunctions H{Z), ipi{Z) and (poiZ) are determined, 
substituting all the fiuctuations (+ background) into the D8-brane action ( |2.9| ) gives 
us the glueball - qq meson couplings. We concentrate on interactions linear in the 
glueball field G, since we are interested in the glueball decays. 

First, because the D8-brane action is written in terms of the type IIA string 
metric and the dilaton field, we need a dimensional reduction from the 11 dimensional 
fields to the 10 dimensional fields. We find 

L f r'^f\ f f F^ \ F^ 

9rr = ;7 + ^^44 + -JjKr j, = — + — /l44 j V^u + -^h^. 



P (l + |^^44 + ^Kr) , e-/3 = ^ + /.44. (2.22) 



7r/t 'Y^i^r^) Qtt jj^. 



■i-Note that there are no higher derivative terms on the right hand side of this expression. This 
is due to a useful gauge choice for the gravity fluctuations (2.4), introduced in | |Tc| ]. 



- 10 - 



Substituting these and all the expressions for the fluctuations (p.4|) and (|2.11|) into 
the D8-brane action (|2.9|) , we obtain the interaction action S'int before the integration 
over the extra dimension Z: 



Tr / d^xdZ 



180K 



(2.23) 



We have defined H{Z) = ((1/4) + 3/{5K - 2))H{Z). In this action (g]^ we have 
kept only terms quadratic in vr and p^, for simplicity. Note that the first line in the 
interaction action (p.23|) vanishes for an on-shell glueball, (□ — M'^)G = 0. 

Terms of higher order in vr and p in the list ( ^l3D can also be computed in the 
same manner. (Other couplings do not appear at this order in the large 't Hooft 
coupling expansion.) Among these additional couplings, only the Gimp coupling 
shown below will be relevant for the later computations of the decay width: 



%s J 



RssU'^iKTr {d,n[p^, vr]) 4>l^,H^G 



(2.24) 



In addition to this, there is another term for Gnnp which vanishes for the on-shell 
glueball. 

Finally, performing the Z integration, we obtain the following interaction la- 
grangian (in this expression, again we have kept only terms quadratic in tt and p^): 



5i 



int 



Tr / rf^x <^ ci-{d^ny 1 



□ 
M2 



G + C2-M^;,p^(l 



^ " " G 



+ci-((9/,7r(9,,7r) 



M2 



-G + C2-Mj^^pf,p^ 



M2 



1 / n \ 1 df^d" 3 (9^ 1 



Here G is the lightest scalar glueball field with the mass M, J^^ = 0~^'^, and F^j^^y = 
dfj,pu — dyPfj, is the field strength of the rho meson p^, and the coefficients ci ~ C4 are 
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defined as follows: 



Cl = / dZ-^H, (2.26) 
J Ktc 

C2 = '^R%'ul/irM{2na'^)g;' J dZ K{dziJifH, (2.27) 

C3 = lR%'u'^irM{2na")g;' J dZ K^'/^^^fH, (2.28) 

C4 = ^<s'f^Ki'^D8V^4(27r«'2)^;iMi^ I dZj^^^Mdz^i)H. (2.29) 



These are evaluated numerically, 



44.3 5.03 



3/2 5 '-'2 ' 

5(ymA^c' Mkk gvuNc' Mkk 

, _ 49.3 _ -0.0732MKK 

C3 — , C4 — jj^ . [A.6\J) 

gyuNc' Mkk 9ymNc' 



The Gnnp coupling (p.24|) is integrated to give 



/d^d^ f 1 ~ 1 43 X 10^ 
rf^x Tr(9^7r[p^, vr])— — G, with 05=/ dZ—^jjiH = — — . (2.31) 

These are the basic ingredients for computing the decay of the lightest glueball in 
section ||. 

2.3.3 Mixing of glueball with qq mesons 

For the identification of the glueball state in the data of the real hadronic spectra, 
mixing with other states possessing the same quantum number (for the scalar glueball 
of our interest, it is 0"*""^) is quite essential. Generically, the mixing is expected to 
appear, because no symmetry can prohibit it. In holographic QCD, mixing can be 
computed explicitly.* Here we show that the lightest scalar glueball has no mixing, in 
the leading order interaction lagrangian of our concern. This means that the mixing 
is largely suppressed, so the decay of the lightest scalar glueball is dominated by a 
direct decay (not through the mixing of the mass matrix). Note that for generic 
glueball excitations, this is not the case, as we will see below. 

First, let us give a generic argument on the mixing of a generic glueball and 
qq mesons, in the holographic QCD. We can show that the mixing is suppressed by 
The mixing, a linear coupling between a glueball and a qq meson, originates 
in a linear coupling between supergravity fields and Yang-Mills/scalar fields in the 
D8-brane action (p.9|). The order of the mixing can be identified after canonically 



* See for glueball mixing to 77' meson in a different holographic model of QCD based on flavor 
D6-branes [ITl. 
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normalizing the glueball field G and a meson field X. We already know that the 
normalization of the glueball field G is given by (|]21|), while that of the X meson 



can be characterized by that of the p meson, ( glSl) . Using the relations (^1201) , the 



pre-factor in (|2.15| ) can be computed as (^ym-^c /108^^- The D8-brane action has the 
same pre-factor (since the pre-factor in ( p.l5|) is basically for canonically normalizing 
the p meson kinetic term in the D8-brane action), so the mixing term can be written 
as 



1 



(fx XG. (2.32) 



This means that the mixing is of order l/y/Wl.'^ 

Although it is suppressed in the large N^. limit, this has a significant effect on 
the decay process. The direct meson decay process comes from meson interactions 
which are of order gyliN'^. So, combining this with the mixing, the total decay 

T 3/2 

amplitude through the mixing is ~ gy-f^Nc . On the other hand, the couplings 
computed in ( p.30| ) mean that the direct decay amphtude of the glueball is of order 



gylj^Nc^^"^, which is the same as the mixing decay amplitude. Therefore in generic 
glueball decay, direct decay process is comparable to the decay through the mixing 
term. 

Our interest here is primarily on the lightest glueball, and let us show that 
there is no mixing for this lightest glueball. First, note that the p meson and the 
pion appear quadratically in the D8-brane action. This already shows that for the 
glueballs originating in the dilaton and the graviton fluctuations have no mixing 
with the p meson and the pion. (For glueballs coming from NSNS 5-field or RR 
gauge fields may have mixings.) So the lightest glueball can mix only with other 
type of mesons which are not in the higher dimensional Yang-Mills field; that is 
the transverse scalar field y on the D8-brane. This y is not written explicitly in 
( p.9| ) but included in the induced metric and the dilaton. The KK decomposition 
of the field y produces scalar mesons with quantum number This y is again 

an Nf X A'^^ matrix. When Nf = 3, among Nj = 9 matrix elements, we have two 
elements with isospin zero, which mix with the glueball. These two qq meson states 
can be identified with /o mesons (other than the glueball candidate /o(1500)); near 
/o(1500), there are /o(1370) and /o(1710), which may be identified with these two 
qq mesons coming from y. 



^This suppression can be understood more easily. The supergravity fields are normalized with 
1 /(/g factor in front of the supergravity action, while the D8-brane gauge fields are normalized with 
1/gs factor in the tension of the D8-brane Tbg. Therefore, if one canonically normalizes the kinetic 
terms of the fluctuation fields, the above-mentioned mixing coupling receives a ^/g^ factor, and in 
view of the AdS/CFT correspondence this factor is just a correction. 
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Possible mixings among these three /o mesons have been studied phenomeno- 
logically (see for example [0, Q). However, we show below that holographic QCD 
predicts there is no mixing at the leading order. -I- From the induced metric in the 
D8-brane action, we are interested in terms linear in the field y, which would lead to 
possible mixing. These are 

9uy\y=od^,y{z,x^'), g^y\y=odzy{z,x^'), gzy\y=od^y{z, xf"), gzy\y=odzyiz,x^), 

y[dygTT,rr,fir,fiu]y=0, y[dy(p]y=0, (2.33) 

where the last term is of course from the dilaton. Here the bulk coordinates y and 



ZUkk are fTl 



rl/2 



f^KK^/^-ljcos^, ^= [t/KK^^-ljsin^, (^^l^^- (2-34) 



It is easy to show that in fact the metrics and the dilaton appearing in ( |2.33|) which 



include the glueball fiuctuations ( |2.4| ) disappear at y = where the D8-brane is 
located, after transforming the metric by using the (r, /i, r) spacetime coordinates. 
So we conclude that there is no mixing of the lightest glueball with mesons, at the 
leading order in Precisely the same argument shows that the lightest 2~^~^ 

glueball also does not participate in mixing at this order. 

3. Decay of the lightest scalar glueball 



Starting from the interaction lagrangian (|2.25| ) and ( p. 311) , we can directly study 



the decay products and their decay widths. In this section, we first enumerate the 
kinematically allowed decay processes by analyzing the masses of particles involved. 
This provides a list of decay products for the lightest glueball. We then compute the 
decay widths by using the interaction lagrangians ( |2.25| ) and ( |2.31| ) including explicit 
numerical coefficients ( |2.30| ). Finally, we compare the widths with experimental 



data for the glueball candidate /o(1500). We find that the prediction of holographic 
QCD qualitatively reproduces the total width as well as the branching ratios of the 
/o(1500). 

3.1 Decay products 

Let us study the kinematical constraints on the decay of the lightest glueball in the 
holographic QCD. The lightest glueball mass M is given by M = ^7.31/9Mkk, 



while the p meson mass is = ^/\^MKK = VO-QGQMkk- So we have a relation 

mp< M < 2mp (3.1) 
^The components with isospin — 1 are identified as ao(1450) meson p7| |. 

^Here we only consider the U{Nf) singlet of the transverse scalars. We can show that the other 
non-Abelian scalars do not mix with the lightest glueball as in the case of the p mesons. 
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in the holographic QCD. This means, our hghtest glueball cannot decay to two on- 
shell p mesons.* 

In ( p.l3| ), we hsted all the coupling terms appearing in the interaction lagrangian 
with a single glueball field G. From those terms one can construct Feynman diagrams 
for the decay processes. We will work with two fiavors, for definiteness. The list ( p.l3| ) 
can be grouped into two categories, as follows: 

(i) GTr(7r2), GTr(7r, [tt, p]), GTr(p2), GH 

(ii) GTr([7r,p]2), GTr(p[p,p]), GTr([p,p]2), Guu 

We have written explicitly and separately the trace part of the qq mesons: rj' for the 
pions and uo for the p mesons. In addition to these, there are couplings coming from 
the Chern- Simons term, 

(iii) G'Tr(7rpp), Gri'Ti^pp), GuTi^np), Grj'uu 

for which the spacetime indices are contracted by the epsilon tensor. The category 
(i) is important for the decay processes, while the category (ii) and (iii) are almost 
irrelevant kinematically, since with the couplings (ii) and (iii) the final decay product 
includes more than 5 pions (or 4 pions and one 77'). To understand this, note that p 
meson can decay to two pions, and u decays to three pions^. When the number of the 
pions are large, typical momentum for the final pion state is small. Pion couplings 
are accompanied by derivatives, then the amplitude is expected to be suppressed. 
By this reason, we restrict our analysis to the cases where the final decay product is 
induced by the couplings (i). 

Since the mass of the glueball is not larger than twice the p meson mass, the 
final decay products should have less than two p mesons. All possible decay chains 
obtained by these couplings are categorized by the decay products: 



(a) 


G- 


TTTT 


(figure 1) 


(b) 


G- 


pmi, 


G —>■ pp —>■ pmr (figure 2) 


(c) 


G- 


pmr - 


— i> TTTTTTTr, G ^ pp ^ TTTTTTTr (figUrC 3) 


(d) 


G- 


-s> 7]'-)]' 


(figure 1) 



If we think of G as the /o(1500), then this list is consistent with what is known in 
the particle data book The branching ratio given in |l| is 35% for (a), 49% for 

*For the most probable candidate of the gluebaU /o(1500), M = 1507MeV, so this mass relation 
is satisfied in the experimental data {nip = 776MeV). 
^We don't consider coupling to photons in this paper. 
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Figure 1: A glueball G decaying to two pions vr. 



G ^ An (corresponding to (b)+(c)-'-), and 7 % for (d).^ So we can reproduce the 
main decay channels of the /o(1500). In the next subsection, we compute the decay 
widths for each of these decay branches. 



3.2 Decay widths 

Let us evaluate the decay widths for these groups. Groups (a) and (d) are two-body 
decays so that the decay widths can be computed analytically. For the remaining 
groups, (b) and (c), integrations over final momenta are complicated that we com- 
puted the decay widths numerically.^ 

3.2.1 G ^ TTTT, G ^ r^T] 

Two-body decays are the simple to analyze, for which we have for the decay width. 



iPl 



57r 



M2 



\M\\ (3.4) 



where M. is the amplitude of the graph responsible for the decay, and p is the final 
momentum of one of the identical particles in the decay product. 



■'■In [Ql, (b) is not explicitly -written, but we interpret that (b) is included in the G — > 47r decay 
in because the on-shell p meson in (b) -would decay to 27r. 

^We are -working in the case of t-wo flavors, so -we don't distinguish and 77', and ignore K . 

^We -worked in the supergravity convention, but to compute the decay -widths it is convenient 
to go to the metric convention (+1,-1,— 1,-1). The ne-w action S — Si^ + Si-at quadratic in the qq 
meson fields is 



^0 = Try d'x |i(a^7r)(a''7r) - \f^.F^^'' + ^AiAfi^p^p^j , 

^int = Try d'x |cii(9^^)(5^^) {l + ^^G + C2 ^Afi^p^p^ (^1 + g 

1 d^d'^ 1 d^^d^ 

-ci-(9^7ra^7r)-^G - C2-M]^j^Pf,p^^^G 



(3.2) 



+cs-F,^F^^ ( 1 _ _ 1 G + c,-F,,F^P^G + c,-p.F^^j^G } . (3.3) 



In the rest frame of the glueball, the first hne in the interaction lagrangian ( p.25| ) 
vanishes. For the 2tt decay, the relevant coupling in that frame is 



^ciaoTT'^aoTr'^G. (3.5) 

For definiteness we consider a specific adjoint index for the pion n"" (a = 1, 2, 3). We 
have two pions as a final state, vr^^-' and 7r^^\ then 

M = ^c^zpf\pf X 2, (3.6) 

where the last factor 2 is for the symmetry of exchanging the two final identical 
particles. The kinematics shows that p^^ = Po = |p| = M/2 because the pions are 
massless, so we obtain 

M = -ciM^A. (3.7) 
The decay width summed over a = 1, 2, 3 is 

- 4^ X 3 X -. (3.8) 

The last factor 1/2 is necessary because the final state has two identical particles. 
For evaluating the numerical value of the decay width, we need k = XNc/ 108n^ = 



7.45 X 10 which was used in |jT^ to fit the pion decay constant, and Nc = 3. Using 



these as inputs, we finally obtain the decay width 



= 0.040. (3.9) 

This is to be compared with the experimental data in [|l|, 

r{^^) 1 no 

= — - X 34.9% = 0.0252, (3.10) 
M 1507 ' ^ ^ 

with which we find a qualitatively good agreement. 

Another two-body decay channel is for G ^ rjrj. The rj' mass evaluated in the 
holographic QCD in |]T^II is found to be too large, 2m^' > M, so in the holographic 



QCD this decay channel cannot be described. However, if we adopt, as a trial, the rj' 
mass as a free parameter in this holographic QCD, then we obtain the decay width 



^G-*riri _ Fg^tttt 1 /-, 



M M 3 \ 



XoX\l-T7T- (3.11) 



'The expression is m^' = 5^\/^(5ym^c)A^kk ~ 17.4AfKK- 
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Figure 2: A glueball G decaying to two pions vr and a single p. There are two graphs, 
the decay with a single vertex (Left) and the decay with two vertices (Right). 



The factor 1/3 is to suppress the effect of the three kinds of the pions, and the last 
factor is necessary to replace |p| of the pion with that of the rj meson. If we substitute 
the real observed ratio ?7i^/iVfj(j(i5oo) = 547.5/1507, we obtain 

= 0.0090. (3.12) 
We compare this with the experimental data. 



(ex) -p(ex) 

/ I I IM 

X 7.0% = 0.00506, (3.13) 



M 1507 
again this is qualitatively in agreement with our result. 

3.2.2 G pmi 

First we describe the decay G prni which uses the single vertex ( p. 311) , see figure 
1^ (Left). The SU{2) generators are normalized as o""/v^, so the interaction ( p.31j ) 
is written explicitly as 

-V2c,ea,,d,7r'^py^G = V2c,eatcdo7r'' p^G (3.14) 
where we have used a relation in the rest frame of the glueball G, 

" " -'^o'^oG'- (3.15) 

Labeling the decay products as 7r-'^(p(i)) + p^(p(p)) + 7r^(p(2)), the amplitude is 

M = V2c5(ip(i)o - «P(2)o)eo (3.16) 

where eo = |p(p)|/mp is the zeroth component of the p meson polarization vector, 

^t, ^ f IP(p)l -P(p)o 
V ' IP(P) 

Other polarization vectors have vanishing eo. 



e'=[^,r^P(p) ■ (3.17) 
V lP(p)l"^p 
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Next, we evaluate the amphtude for the process G ^ pp ^ pim, see figure || 
(Right). Let us hst the Feynman rules: The n'^p^Tr'^G vertex is the same as before, 
V^csip'^ ■'^o^afec, and the vr'^p^vr-^ vertex was obtained in as 

r°° 1 24 

\Pic^iV^f'>^def where Cg = / dZ—-^pi = —— — . (3.18) 



oo 



The p'tpIG vertex (no sum over a) is 



and the p meson propagator is 



Pfp) - "^p + 'in^'pTp y m 

Here we need the p meson decay width, Tp/nip = Cg/247r, which can be evaluated as 
Tp/nip = 0.307 in the holographic QCD [17, 18|. This is close to the experimental 
value Tp/mp = 149.4/775.5 = 0.1927 in |1|]. As for the polarization vectors, in 
addition to ( |3.17| ), we have two more vectors. One is 

e^= \0, I (3.21) 



\/ip['^r + ipi'^y \/ip['^r + ipi'^y 



and the other gives the same decay width with that of the polarization (|3.21| ). Using 
all of these Feynman rules, we can compute the decay amplitude for the process 
G — s> pp — s> pun. 

The total expression for the decay width is lengthy, and we here provide only the 
numerical results after substituting the necessary inputs used also in the evaluation 
of the Fc^TrTr- We obtain 



G — >p7Tn 

M 



3 X (l.l X 10"^ + 2 X (2.1 X 10"^)) ~ 1.3 x 10"^ (3.22) 



The factor 3 is for the sum over possible fiavor indices. (The decay width we com- 
puted is for the decay G Tr^p^vr^, and there are two other possible combinations 
for the indices.) Terms in the parentheses are for different polarizations of the final 
p meson. 

Our result (|3.22|) is very small. The smallness mainly comes from the fact that 
the integration region of the momentum is so small because rup is very close to M 
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Figure 3: A glueball G decaying to four pions tt. There are two graphs, the 
decay with two vertices (Left) and the decay with three vertices (Right). 



in our holographic computation. In reahty, the mass of /o(1500) is much larger than 
the mass of the p meson. So, as a trial, in our computation of the decay width, let us 
modify the input glueball mass M such that M/nip coincides with the experimental 
value in Then, we obtain 

^ G—tpTTTT 



M 



0.096. 



(3.23) 



Since the decay product pvrvr seems to be included in G — 47r in [|I| , we compare our 
result with the experimental value after adding the decay width oi G ^ Att which 
we compute next. 

3.2.3 G ^ 47r 

The computation of this amplitude is done in the same manner, and we do not write 
it explicitly, except for some important points. First, it is easy to find out that the 
decay product is only in the combination G 27r*27r-^ {i ^ j). So one can specifically 
choose i = 1, j = 2 for the computation. (This determines the index for the p meson 
as p^.) Since the amplitude is proportional to 



all possible ways to assign the adjoint index for each of the final pions are 



(3.24) 



7r^(P{i)) 


a 


d 


a 


d 


a 


f 


a 


f 


c 


f 


c 


f 


c 


d 


c 


d 


^^(P(2)) 


d 


a 


d 


a 


f 


a 


f 


a 


f 


c 


f 


c 


d 


c 


d 


c 


^'(P{3)) 


c 


c 


f 


f 


c 


c 


d 


d 


a 


a 


d 


d 


a 


a 


f 


f 


7r^(P(4)) 


f 


f 


c 


c 


d 


d 


c 


c 


d 


d 


a 


a 


f 


f 


a 


a 


sign 


+ 


+ 


+ 


+ 










+ 


+ 


+ 


+ 











We have to sum the amplitude with these substitutions, with the sign indicated in 
the table. 

The main difficulty in the computation resides in the evaluation of the integration 
of the final momenta, / d?^{i)d?Y>{2)d^V{?,)d^V{i) ■ the integrand of the decay width. 
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there is a four dimensional delta function coming from the total energy momentum 
conservation. The integration J (i^p(4) trivially eliminates three delta functions of 
the momentum conservation. Furthermore, using spatial rotation symmetry, we 
can orient one of the remaining momenta as p(3) = (p, 0,0), then the remaining 
single delta function for the energy conservation can be eliminated by the integration 
/ d^P(3) = 47r / p^dp. Specifically, the delta function is expressed as 



S {P+ \Ip^ + MP(1)x + P{2)x) + |P(1) + P(2)P + |P(1)| + |P(2)| 



M 



2(M- 


P(i) 


- |P(2)|)(P(l)x 


+ P(2)x) 


+ 1 


P(l) + P(2) 


2 + (M- 


P(l) 




P(l)l)^ 


2(M- 


P(l) 




P(2)| + P{l)x + P{2)xy 



(M-|] 




- |] 


^(2)1)^- 


P(l)+P(2)|' 


2(M- 


|P(i 




|P{2) 


+ P{l)x + P{2)x) 



x6[p - ' '7^^' , '7^^'^, . (3.25) 



So the integration over p results in a constraint for the remaining momenta P(i) 
and P(2). This constraint corresponds to a restriction on the integration region, 

|P(1)I + |P(2)| + |P(1) + P(2)| < M. 

Numerical integration for the remaining momenta gives the decay width,* 

- 2.2 X 10-^ (3.26) 



M 



If we adjust the glueball mass to the experimental value (while fixing the p meson 
mass by nip = y/XiM^K) as was done before, we obtain 



r 



~ 0.0087. (3.27) 



M 



Let us compare our results (|3.23| ) and ( p.27| ) with the experimental values. 
Adding these two, we obtain 

^G^4n + TG^p^^^_^^^^ (3.28) 



M 

while the experimental data shows 

= X 49-5% = 0.0358. (3.29) 

One can see that the order of this decay width is reproduced in the holographic 
QCD.t 



*This value includes a factor 3 accounting for different combinations of the species of the final 
decay product. 



Hi we include the mass of the pions in some way, the decay width (3.28) to the four pions is 



expected to become significantly smaller and close to the experimental data. 
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4. Concluding remarks 



We have presented here the first attempt in computing decays of glueballs to qq 
mesons using holographic QCD. We have adopted a string-theoretic set-up, (which 
is of the so-called "top-down" type), the Sakai-Sugimoto model. The glueball sector 
lives in supergravity fluctuations in the Witten's background of non-BPS black 4- 
branes, and the mesons live on the probe D8-branes. The coupling between the two 
sectors is encoded in the D8-brane action, and KK decomposition and integration 
over extra dimensions gives the desired couplings in four spacetime dimensions. 

Explicit couplings between the lightest glueball and the qq mesons are given, and 
the associated decay products/widths are calculated. We find that our results are 
consistent with the experimental data of the decay for the /o(1500) which is thought 
to be the best candidate of a glueball in the hadronic spectrum. 

The most important merit of the holographic QCD is that one can go beyond 
the chiral perturbation theory; one can compute coefficients which cannot be fixed 
solely by the chiral symmetry. At low energy the chiral perturbation works well, 
but at the energy scale of the glueball mass the derivative expansion in the chiral 
perturbation becomes unreliable. Furthermore, glueballs are fiavor-blind, so it is 
quite difficult to constrain possible interactions from the chiral symmetry. The holo- 
graphic description obtained in the holographic QCD is, in principle, equivalent to 
QCD, though in the large and large 't Hooft coupling limit. We therefore ex- 
pect that the holographic approach should provide interesting information on strong 
coupling physics of QCD. In fact, we have discussed generic features of the glueball 
interactions predicted in holographic QCD (see section |2.2|) . For example, we have 
argued that, among the decay products of glueballs, Atiq should be suppressed. 

One of the reasons why the /o(1500) is expected to be a glueball state is that the 
/o(1500) does not decay to 27. In the holographic QCD, we can compute relevant 
photon coupling in the same manner, and find that G77 coupling is vanishing at 



the leading order (see |jT^ ^ for the way to introduce the electromagnetic field 
as an external background of the massless fields on the D8-branes). Since we have 
shown that there is no mixing with qq mesons at the leading order, our result of the 
suppressed photon coupling reproduces the experimental data. 

In this paper, we have explicitly computed for the decay of the lightest glueball, 
which is of the most phenomenological interest. There are also many other interest- 
ing directions, e.g., generalizing our results using various approaches in holographic 
QCD. Here are some examples: 



■i-Precisely speaking, the limit to QCD includes A/kk ^ 00 (there should be a double-scaling 
limit with a simultaneous scaling of the 't Hooft coupling), after incorporating infinite number 
of 1/Nc corrections. The background receives large stringy corrections and becomes essentially a 
purely stringy background. 



- 22 - 



Multi-glueball couplings. Self-couplings of the glueballs can be computed in the 
supergravity sector. Emission of mesons from a propagating glueball can be 
described by the D8-brane action similarly. For highly spinning glueballs whose 
holographic dual are closed fundamental strings in the confining supergravity 



background, their decay into two glueballs was briefly described in [22 



Universally narrow width of glueballs. If one can show in the holographic 
QCD that the total decay width of any glueball state is narrow, that would 
provide support for this widely-held belief. In this paper we have shown the 
narrowness only for the lightest glueball. Explicit calculation of the widths is 



possible for other glueball excitations, as they are available in |jT^. The 2"*"+ 
glueball coupling has been described in (|2.14| ) for example. 

• Glueballs with other J"^*" quantum numbers. Glueball states originating in the 
RR fields in the supergravity may possess interesting structure in the meson 
couplings. The 0~+ glueball is described by a RR 1-form whose fluctuation 
is completely decoupled from the others, and it appears in the Chern-Simons 
coupling in the D8-brane action. The 1^ glueballs reside in the NS-NS 2- 
form field, and it should have a large mixing with the meson fields. This is a 
consequence of the gauge invariance in the supergravity, requiring the gauge- 
invariant combination -Bnsns + F m the D8-brane action. 

• Thermal/dense QCD. One can modify the supergravity background, or in- 
troduce a background for the D8-brane fields, to describe the finite temper- 
ature/baryon density^, that surely will modify the couplings which we have 
computed. Glueball couplings should be sensitive to the deconfinement tem- 
pearture, and, near the transition temperature, they should become singular 
in some sense. This is of phenomenological interest in view of the onset of the 
LHC. 

• Computation of the glueball couplings in other models of holographic QCD. For 
a single-flavor case, flavor D6-branes enable one to introduce easily the quark 
mass [l^, which might shed some light on how our results may be modified 
by the pion masses. To apply our strategy to so-called bottom-up phenomeno- 
logical approach in holographic QCD may reveal how universal the glueball 
couplings obtained in our paper are. 

All of these would be interesting to investigate. 



§ There are many recent works on introducing finite temperature or finite baryon density (or 
chemical potential) in holographic models of QCD. Specifically, for the Sakai-Sugimoto model, see 



- 23 - 



Acknowledgments 



K. H. would like to thank T. Sakai and S. Sugimoto for valuable comments, and like to 
thank Brown university high energy theory group for kind invitation. S. T. is grateful 
to D. Jido and M. Murata for useful discussions. K. H. and S. T. thank the Yukawa 
Institute for Theoretical Physics at Kyoto University which supported the YITP 
workshop YITP-W-07-05 on "String Theory and Quantum Field Theory" during 
which valuable discussions on the present paper were made. K. H. and S. T. are 
partly supported by the Japan Ministry of Education, Culture, Sports, Science and 
Technology. 

References 

[1] W. M. Yao et al. [Particle Data Group], "Review of particle physics," J. Phys. G 33, 
1 (2006). 

[2] C. Amsler and N. A. Tornqvist, "Mesons beyond the naive quark model, " Phys. 
Kept. 389, 61 (2004). 

[3] C. J. Morningstar and M. J. Peardon, "The Glueball spectrum from an anisotropic 
lattice study," Phys. Rev. D 60, 034509 (1999) [arXiv:hep-lat/9901004]. 

C. Michael, "Hadron spectroscopy from the lattice, " AIP Conf. Proc. 432, 657 
(1998) [arXiv:hep-ph/9710502]. 

[4] W. J. Lee and D. Weingarten, "Scalar quarkonium masses and mixing with the 
lightest scalar glueball," Phys. Rev. D 61, 014015 (2000) [arXiv:hep-lat/9910008]. 

[5] J. M. Maldacena, "The large N limit of superconformal field theories and 

supergravity," Adv. Thcor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 
(1999)] [arXiv:hep-th/9711200]. 

[6] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "Gauge theory correlators from 
non-critical string theory, " Phys. Lett. B 428, 105 (1998) 
[arXiv:hep-th/9802109]. 

E. Witten, "Anti-de Sitter space and holography, " Adv. Theor. Math. Phys. 2, 253 
(1998) [arXiv:hep-th/9802150]. 

[7] E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge 
theories," Adv. Theor. Math. Phys. 2, 505 (1998) [arXiv:hep-th/9803131]. 

[8] C. Csaki, H. Ooguri, Y. Oz and J. Terning, "Glueball mass spectrum from 
supergravity," JHEP 9901, 017 (1999) [arXiv:hep-th/9806021]. 

R. de Mello Koch, A. Jevicki, M. Mihailescu and J. P. Nunes, "Evaluation of 
glueball masses from supergravity," Phys. Rev. D 58, 105009 (1998) 
[arXiv:hep-th/9806125]. 



- 24 - 



[9] D. J. Gross and H. Ooguri, "Aspects of large N gauge theory dynamics as seen by 
string theory," Phys. Rev. D 58, 106002 (1998) [arXiv:hep-th/9805129]. 

M. Zyskin, "A note on the glueball mass spectrum," Phys. Lett. B 439, 373 (1998) 
[arXiv:hep-th/9806128]. 

H. Ooguri, H. Robins and J. Tannenhauser, "Glueballs and their Kaluza-Klein 
cousins," Phys. Lett. B 437, 77 (1998) [arXiv:hep-th/9806171]. 

J. G. Russo, "New compactifications of supergravities and large N QCD, " Nucl. 
Phys. B 543, 183 (1999) [arXiv:hep-th/9808117]. 

A. Hashimoto and Y. Oz, "Aspects of QCD dynamics from string theory," Nucl. 
Phys. B 548, 167 (1999) [arXiv:hep-th/9809106]. 

C. Csaki, Y. Oz, J. Russo and J. Terning, "Large N QCD from rotating branes," 
Phys. Rev. D 59, 065012 (1999) [arXiv:hep-th/9810186]. 

J. A. Minahan, "Glueball mass spectra and other issues for supergravity duals of 
QCD models," JHEP 9901, 020 (1999) [arXiv:hep-th/9811156]. 

J. G. Russo and K. Sfetsos, "Rotating D3 branes and QCD in three dimensions, " 
Adv. Theor. Math. Phys. 3, 131 (1999) [arXiv:hep-th/9901056]. 

C. Csaki, J. Russo, K. Sfetsos and J. Terning, "Supergravity models for 3+1 

dimensional QCD," Phys. Rev. D 60, 044001 (1999) [arXiv:hep-th/9902067]. 

[10] N. R. Constable and R. C. Myers, "Spin-two glueballs, positive energy theorems and 
the AdS/CFT correspondence," JHEP 9910, 037 (1999) [arXiv:hep-th/9908175]. 

[11] R. C. Brewer, S. D. Mathur and C. I. Tan, "Discrete spectrum of the graviton in the 
AdS(5) black hole background," Nucl. Phys. B 574, 219 (2000) 
[arXiv:hep-th/9908196]. 

[12] R. C. Brower, S. D. Mathur and C. I. Tan, "Glueball spectrum for QCD from AdS 
supergravity duality," Nucl. Phys. B 587, 249 (2000) [arXiv:hep-th/0003115]. 

[13] A. Karch and E. Katz, "Adding flavor to AdS/CFT, " JHEP 0206, 043 (2002) 
[arXiv:hep-th/0205236]. 

[14] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, "Meson spectroscopy in 
AdS/CFT with flavour," JHEP 0307, 049 (2003) [arXiv:hep-th/0304032]. 

T. Sakai and J. Sonnenschein, "Probing flavored mesons of confining gauge theories 
by supergravity," JHEP 0309, 047 (2003) [arXiv:hep-th/0305049]. 

J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik and L Kirsch, "Chiral 

symmetry breaking and pions in non- super symmetric gauge / gravity duals, " 
Phys. Rev. D 69, 066007 (2004) [arXiv:hep-th/0306018]. 

P. Ouyang, "Holomorphic D7-branes and flavored N = 1 gauge theories," Nucl. 
Phys. B 699, 207 (2004) [arXiv:hep-th/ 03 11084]. 



- 25 - 



C. Nunez, A. Paredes and A. V. Ramallo, "Flavoring the gravity dual of N = 1 
Yang-Mills with probes," JHEP 0312, 024 (2003) [arXiv:hep-th/0311201]. 

S. Hong, S. Yoon and M. J. Strassler, "Quarkonium from the fifth dimension," 
JHEP 0404, 046 (2004) [arXiv:hep-th/0312071]. 

N. J. Evans and J. P. Shock, "Chiral dynamics from AdS space," Phys. Rev. D 70, 
046002 (2004) [arXiv:hep-th/0403279]. 

B. A. Burrington, J. T. Liu, L. A. Pando Zayas and D. Vaman, "Holographic duals 
of flavored N = 1 super Yang-Mills: Beyond the probe approximation," JHEP 
0502, 022 (2005) [arXiv:hep-th/0406207]. 

R. Casero, C. Nunez and A. Paredes, "Towards the string dual of N = 1 SQCD-like 
theories," Phys. Rev. D 73, 086005 (2006) [arXiv:hep-th/0602027]. 

[15] D. T. Son and M. A. Stephanov, Phys. Rev. D 69, 065020 (2004) 
[arXiv:hep-ph/0304182]. 

J. Erhch, E. Katz, D. T. Son and M. A. Stephanov, "QCD and a holographic model 
ofhadrons," Phys. Rev. Lett. 95, 261602 (2005) [arXiv:hep-ph/0501128]. 

L. Da Rold and A. Pomarol, "Chiral symmetry breaking from five dimensional 
spaces," Nucl. Phys. B 721, 79 (2005) [arXiv:hep-ph/0501218]. 

A. Karch, E. Katz, D. T. Son and M. A. Stephanov, "Linear confinement and 
AdS/QCD," Phys. Rev. D 74, 015005 (2006) [arXiv:hep-ph/0602229]. 

G. F. de Teramond and S. J. Brodsky, "The hadronic spectrum of a holographic 
dual of QCD," Phys. Rev. Lett. 94, 201601 (2005) [arXiv:hep-th/0501022]. 

K. Ghoroku, N. Maru, M. Tachibana and M. Yahiro, "Holographic model for 
hadrons in deformed AdS(5) background," Phys. Lett. B 633, 602 (2006) 
[arXiv:hep-ph/0510334]. 

[16] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, "Towards a holographic 
dual oflarge-N(c) QCD," m¥P 0405, 041 (2004) [arXiv:hep-th/0311270]. 

[17] T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD, " Prog. 
Theor. Phys. 113, 843 (2005) [arXiv:hep-th/0412141]. 

[18] T. Sakai and S. Sugimoto, "More on a holographic dual of QCD," Prog. Theor. 
Phys. 114, 1083 (2006) [arXiv:hep-th/0507073]. 

[19] K. Hashimoto, T. Hirayama and A. Miwa, "Holographic QCD and pion mass," 
JHEP 0706, 020 (2007) [arXiv:hep-th/0703024]. 

O. Bergman, S. Scki and J. Sonnenschein, "Quark mass and condensate in HQCD," 
arXiv: 0708 .2839 [hcp-th] . 

A. Dhar and P. Nag, "Sakai- Sugimoto model, Tachyon Condensation and Chiral 
symmetry Breaking," arXiv:0708.3233 [hep-th]. 



- 26 - 



[20] J. L. F. Barbon, C. Hoyos, D. Mateos and R. C. Myers, "The holographic life of the 
eta'," JHEP 0410, 029 (2004) [arXiv:hep-th/0404260]. 

[21] C. Amsler and F. E. Close, "Evidence for a a scalar gluehall," Phys. Lett. B 353, 
385 (1995) [arXiv:hep-ph/9505219]. 

D. Weingarten, "Scalar quarkonium and the scalar gluehall," Nucl. Phys. Proc. 
Suppl. 53, 232 (1997) [arXiv:hep-lat/9608070]. 

F. E. Close and A. Kirk, "Scalar gluehall q anti-q mixing ahove 1-GeV and 
implications for lattice QCD," Eur. Phys. J. C 21, 531 (2001) 
[arXiv:hep-ph/0103173]. 

F. Giacosa, T. Gutschc, V. E. Lyubovitskij and A. Faessler, "Scalar meson and 
gluehall decays within a, effective chiral approach, " Phys. Lett. B 622, 277 
(2005) [arXiv:hcp-ph/0504033]; "Scalar nonet quarkonia and the scalar gluehall: 
Mixing and decays in an effective chiral approach, " Phys. Rev. D 72, 094006 
(2005) [arXiv:hep-ph/0509247]; 

A. H. Fariborz, "Isosinglet scalar mesons helow 2-GeV and the scalar gluehall 

mass," Int. J. Mod. Phys. A 19, 2095 (2004) [arXiv:hep-ph/0302133]; "Prohing 
the Quark and Gluehall Admixtures of the Scalar Mesons, " Int. J. Mod. Phys. 

A 19, 5014 (2004); "Mass uncertainties of f0(600) and f 0(1370) and their 
effects on determination of the quark and gluehall admixtures of the 1=0 
scalar mesons," Phys. Rev. D 74, 054030 (2006) [arXiv:hep-ph/0607105]. 

[22] A. L. Cotrone, L. Martucci and W. Troost, "String splitting a,nd strong coupling 
meson decay," Phys. Rev. Lett. 96 (2006) 141601 [arXiv:hep-th/0511045]. 

[23] O. Aharony, J. Sonnenschein and S. Yankielowicz, "A holographic model of 
deconfinement and chiral symmetry restoration," Annals Phys. 322, 1420 
(2007) [arXiv:hep-th/0604161]. 

N. Horigome and Y. Tanii, "Holographic chiral phase transition with chemical 
potential," JHEP 0701, 072 (2007) [arXiv:hep-th/0608198]. 

K. Y. Kim, S. J. Sin and I. Zahed, "The Chiral Model of Sakai-Sugimoto at Finite 
Baryon Density," arXiv:0708.1469 [hep-th]. 



-27- 



